/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phases (c4 water);

radiationModel none;

c4
{
    type detonating;
    // values for c-4 from: Alia2006
    reactants
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState Murnaghan;
        }
        equationOfState
        {
            Gamma           0.25;
            pRef            1e5;
            rho0            1770;
            K0              0;
        }
        specie
        {
            molWeight       28.97;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1400;           // Heat capacity
            Hf              0.0;
        }
    }
    products
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState JWL;
        }
        equationOfState
        {
            rho0      1770;
            A         6.17e11;
            B         0.16926e11;
            R1        4.4;
            R2        1.2;
            omega     0.25;
        }
        specie
        {
            molWeight       55.0;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1400;           // Heat capacity
            Hf              0.0;
        }
    }

    activationModel linear;
    initiation
    {
        E0          10.1e9;
        points     ((0 0.1143 0));
        vDet       8300;
    }

    residualRho     1e-6;           // Minimum density of the phase
    residualAlpha   1e-10;          // Minimum volume fraction used for division
}

water
{
    type basic;
    thermoType
    {
        transport   const;
        thermo      eConst;
        equationOfState linearTillotson;
    }
    equationOfState
    {
        p0      1e5;
        rho0    1000.0;
        e0      3.542e5;
        omega   0.28;
        A       2.2e9;
        B       9.94e9;
        C       14.57e9;
        pCav    5000;
    }
    specie
    {
        molWeight       18;
    }
    transport
    {
        mu              0;
        Pr              1;
    }
    thermodynamics
    {
        type        eConst;
        Cv          2400; // Give correct reference temperature
        Hf          0;
    }

    residualRho     1e-6;
    residualAlpha   1e-10;
}

// ************************************************************************* //
